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ABSTRACT 

We present a study of the properties of the 0-type, massive eclipsing binary 2MASS J13130841-6239275 located in the outskirts of 
the Danks 2 cluster in the G305 star-forming complex, using near-infrared spectroscopy from VLT/ISAAC. We derive the masses 
and radii to be 24.5 ± 0.9 Mo and 9.2 ± 0.1 R© for the primary and 21.7 ± 0.8 Mq and 8.7 ± 0.1 R© for the secondary component. In 
addition, we evaluate the sensitivity of our parameters to the choice of the spectral features used to determine the radial velocities. Both 
components appear to be main-sequence 06.5-07 type stars at an age of ~ 5 Myr, which is in agreement with the age of the cluster. 
A high visual extinction of A5495 = 11.9 ± 0.1 mag is reported, which is likely attributed to the cold molecular gas contaminating the 
north-east region of the cluster. By fitting the spectral energy distribution of the system to the available BVIcJHK, photometry, we 
determine a distance to the system of 3.52 ± 0.08 kpc with a precision of 2%, which is the most well-determined distance to the Danks 
2 cluster and the host complex reported in the literature. 

Key words, binaries: eclipsing - open clusters and associations: individual: Danks - stars: distances - stars: early-type - stars: 
massive 


1. Introduction 


The G305.4H-0.1 complex (hereafter, G305) is one of the 
most massive, ongoin g star-forming regions in the Galaxy 

HII regions 


Clark &Port9 |2004. It s home to nine _ _ 

Caswell & HavneJ 1987tl . which are powered by the ionizing 
flux of hot stars and are spatially associated with m aser emis¬ 
sion and at least five star clusters dPutra et al.ll2003h embedded 
in the complex. Whereas the majority of these clust ers are com¬ 
pact and thu s unresolved, the Danks 1 & 2 clusters (iDanks et al.l 
19831 Il984ll . which reside in the center of the complex, are re¬ 
solved and thus allow for studies of their physical properties. 
Both clusters have been shown to be < 7 Myr old, with D anks 2 
being noticeably older than Danks 1 dChene et al.ll2012h . In ad¬ 
dition, initial masses of 8 000 and 3 000 Mq were derived for 
Danks 1 and 2 respectivel y, with an uncertainty that reaches 
~ 30% dPavies et al.ll2012l) . Constraining the properties of the 
clusters relies on an accurate determination of the distance, ex¬ 
tinction and their associated uncertainties. To date, the most ro¬ 
bust method to provide accurate distance measurements to young 
galactic clusters an d nearby galaxies is wit h the use of eclipsing 
binaries (EBs) fe.g. lPietrzvns^ki et aDl2013l) . 


* Based on observations made with ESO Telescopes at the La Silla 
Paranal Observatory under program ID 090.D-0065(A). 

** mkourniotis@astro.noa.gr 


EBs exhibiting double-lined spectra provide a geometrical 
method to pr ecisely measure the fundamental para meters of their 
components dAndersenlll99ll : [Torres et ^120101) . In particular, 
the light curve provides the orbital period, eccentricity, the flux 
ratio of the two stars, fractional radii and the inclination of the 
system with respect to the observer. Erom the double-lined spec¬ 
tra, we measure the velocities of the individual components. The 
fit to these velocities yields the velocity semi-amplitudes, and 
therefore, the ratio of the masses. Effective temperatures and thus 
luminosities can be estimated by fitting synthetic spectra to the 
observed ones. Therefore, an accurate determina tion of the red¬ 
dening law and thereby distance, is feasible (e.g. lBonanos et all 
I 2 OO 6 I I 20 II]) . Selecting systems that host well-separated com¬ 
ponents (i.e. detached systems), avoids complexities that orig¬ 
inate from a common-envelope state. Such systems are of ma¬ 
jor importance for testing theoretical stellar evolutionary models 
of single stars. Light curves of detached systems that display 
eclipses of similar depth ensure roughly equal contribution of 
both components to the total flux and such systems with orbital 
periods <10 days are likely to provide double-lined and well- 
separated spectra. 

The majority of hot, massive stars are claimed to be in close 
binar y systems dChini et al.l 1201 2t ISana et aP 120121 : ISota et all 
I 2 OI 4 I) . thus increasing the po ssibility of observing s hort-period, 
eclipsing configurations (e.g. iKourniotis et al.ll20l4) . Such sys¬ 
tems not only constitute extragalactic anchors, ideal for the cal- 
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ibration of th e cosm i c distance scale (e . g. iGuinan et al J 1998 
Hilditch et all l200l [Rib as et a I] 120051; iBonanos et al.l 1200^ 
20 lH) . but also, when found in young, Galactic clusters, can 
set tight co nstrains on the age, chemical composition and dis¬ 
tance (e.g. ISouthworth & ClausenI I20d5 iKoumpia & BonanosI 
I 2 OI 2 I) . In this work, we investigate the properties of 2MASS 
J13130841-6239275 (hereafter, D2-EB), a massive EB located 
~2' from the center of the Banks 2 cluster that constitutes a lu¬ 
minous, double-lined spectroscopic system. A summary of the 
parameters of D2-EB derived in the present work along with the 
parameters of the Banks 2 cluster, are listed in Table [1] 

B2-EB0] was discovered to be a well-detached eclipsing sys¬ 
tem during the optical variability survey by Bonanos et al. (in 
prep.). Bue to the heavy foreground visua l extinction in th e di¬ 
rection of the clusters (Ay = 8-9 mag. iBica et al.l 1^0041) . we 
obtained near-infrared /T-band spectroscopy to resolve double 
lines in the system. The paper is organized as follows: in Sect. |2] 
we describe the follow-up spectroscopy and the data reduction, 
in Sect. [3 we present the radial velocity analysis, in Sect.|4]we 
model the binary, in Sect. |5]we derive the distance, in Sect. |6]we 
discuss our results and in Sect.|2]we finish with some concluding 
remarks. 


2. Observations and data reduction 

Bonanos et al. (in prep.) undertook an optical variability survey 
of the Banks clusters in search of massive EBs over 25 nights in 
March, April and May 2011 with the 1 -m Swope telescope at Las 
Campanas Observatory in Chile. The telescope is equipped with 
a 2048 X 3150 detector with a pixel scale of 0.435"/pixel and a 
field of view of 14.8' x 22.8'. Observations were taken between 
HJB 2455 648.5 and 2455 687.83, in the broadband B,V,R,Ic 
filters over approximately 50, 500, 950 and 530 epochs, respec¬ 
tively. The data were reduced, calibrated to the standard system 
and the ima ge subtraction package ISIS (lAlard & Luptonlll998l: 

ardll2000l) was used to extract the /c-band light curves of all 
sources. The light curves were searched for periodicity with the 
Analy sis of Variance (AoV) algorithm (ISchwarzenberg-Czemvl 
Il989ll yielding 21 EBs. Of these, 13 were found to display red 
colors of R - Ic > 0, likely corresponding to members of the 
highly extincted clusters. The /c-band light curve of B2-EB was 
found to display a strong periodic signal of ~ 3.37 days and 
appeared to be an ideal target for follow up spectroscopy: a lu¬ 
minous, well-detached and nearly equal-depth eclipsing system. 

2.1. Spectroscopy 

We obtained ^f-band spectroscopy of the B2-EB with the In¬ 
frared Spectrometer and Array Camera (ISAAC) mounted at the 
Nasmyth A focus of the 8.2m UT3 (Melipal) at the VET facil¬ 
ity in Chile. The observations were taken under ESO program 
IB 090.B-0065(A). The spectrograph is equipped with a 1024 x 
1024 Hawaii Rockwell array with a pixel scale of 0.148"/pixel 
(iMoorwood et alJ[T998l) . A slit width 0.3" was used, resulting in 
a resolving power of R ~ 7 500 as measured from the full width 
half maximum (~ 2.88 A) of an unblended telluric OH emis¬ 
sion feature at ~ 2.125 pm. To derive a precise measurement 
of the velocity semi-amplitudes, we requested eight visits, four 
per quadrature. In total, eight observing blocks were executed in 
service mode, six of which correspond to the first and two to the 
second quadrature of the system. Eor each observing block, nod- 


‘ RA = 13:13:08.41, Dec = -62:39:27.5, J2000 


ding 60" along the slit was applied in an AB pattern, resulting in 
two spectra of total integration time 100 sec each. 

The spectra were reduced using basic IRAE routine^ 
and wavelength cal ibrated using telluric OH sky lines 
(iRousselot et al.ll2000l) that span the fiT-band wavelength range. 
Subsequently, the extracted pair of spectra of each observing 
block were normalized and averaged, yielding a signal-to-noise 
ratio (S/N) ~ 150 as measured in the middle of the wavelength 
coverage (~ 2.14 ^um). Cosmic rays were rejected using the 
L.A.Cosmic script (Ivan Bok kuml200 i]). We employed the newly 
released molecfit tool dSmette et al.T2015h . which generates a 
synthetic model of atmospheric transmission based on the infor¬ 
mation provided by the EITS header of the science object, to re¬ 
move telluric features. Heliocentric corrections were calculated 
using IRAE’s rvcorr task and applied. Table |3 lists the Helio¬ 
centric Julian Bates (HJB) of the observations, exposure times 
and the values of airmass. Of the eight spectra, we discarded the 
one obtained on HJB 2456 378.80874 from further analysis, as 
it was taken near an eclipse, at phase ~ 0.47, and contained only 
one set of spectral features. 


3. Radial velocity analysis 

The spectrum of the B2-EB clearly displays He ii (2.189 pm) in 
absorption, which is prominent for both components in phases 
outside of eclipse. As in the optical, Heii lines appear only in 
hotter stars and are the main indicator for identifying O-type 
stars. The observed Hei (2.112 yum) in absorption implies that 
both stars must be later than 05 (iHanson et al.lll99^ . Both com¬ 
ponents of the B2-EB are hence classified as intermediate/late 
O-type stars. In addition, our spectra show Bry (2.166 pm) in 
absorption, which provides evidence against a supergiant clas¬ 
sification. However, a precise luminosity classification strongly 
depends on the resolution and requires a S/N > 150. Emission 
of C IV (2.08 pm) is evident only for the redshifted component 
as the blueshifted feature is not within the observed wavelength 
range. Eurther analysis of the observed features that are used to 
measure the radial velocities and to assign a spectral type is pre¬ 
sented in Sec. 13.21 

3.1. Analysis Using Model Atmospheres 

We used the stellar a tmosphere code EASTWINB 
(ISantolava-Rev et all 119971: iPuls et al] l2005l) to synthesize 
the lines of H and He and measure the radial velocities. The 
code enables non-LTE calculations and assumes a spherical 
symmetry geometry. The wind velocity structure is included in 
the model through a /J-like law. We generated a grid of EAST¬ 
WINB templates at the metallicity of the Galaxy, spanning the 
temperature range 35^0 kK with a step of 1 kK and fixed the 
surface gravity at log g - 3.9 dex, as inferred from Sect.|4] The 
templates were broadened and resampled to match the resolution 
{R ~ 7 500) and dispersion of the observed spectra. To constrain 
the projected rotational broadening, we implemented a code 
based on minimization setting six free parameters, three for 
each component; temperature, radial and projected rotational 
velocity. We built synthetic spectra and fit the regions around 
the Bry (2.166 pm)/ Hei (2.162 pm) blend and Hen (2.189;um) 
feature of the ISAAC spectrum taken on HJB 2456 345.85671 

^ IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 
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that exhibits well-separated spectral features. Each composite 
spectrum consists of two FASTWIND templates weighted 
according to the light ratio provided by the light curve analysis 
and broadened at the projected rotational velocity range 50 - 200 
km s '. For the fit, we adopted a cr value of 0.0067, which 
corresponds to the standard deviation of the noise in units of 
normalized flux, based on the mean S/N ratio ~ 150. Fig. [T] 
displays a panel of our analysis providing the reduced-;^f^ of 
the fit with respect to the temperatures and projected rotational 
velocities of the components. For each set of parameters, radial 
velocities were measured to optimize the fit. We found that the 
projected rotational velocities that correspond to the best fit are 
150 + 20 km s ' and 130 + 20 km s ' for the primary and the 
less luminous secondary component, respectively. These values 
are in good agreement with those of a tidally locked system 
whose orbital period is identical to the rotational period for both 
components (Sect.|4]i. 

To construct the radial velocity curve, we fixed the afore¬ 
mentioned projected rotational velocities and ran the fitting algo¬ 
rithm for the seven ISAAC spectra. Given that the best-fit tem¬ 
peratures are found to range from 36 kK to 38 kK with a 1 kK 
difference between the temperature of the components, we even¬ 
tually examined two cases: a set with 37 kK and 38 kK models 
and a set with 36 kK and 37 kK models. While the former re¬ 
sulted in a best fit for the He ii line, the latter appeared to opti¬ 
mally fit the Bry/He i feature. We adopted the first case for pro¬ 
viding reliable velocity measurements as He ii is a more accurate 
tracer of the stellar photosphere for stars at these temperatures. 
Contours of values for the He ii fit were plotted in two di¬ 
mensional radial velocity diagrams and the center of the contour 
corresponding to the lowest was adopted to generate the best 
fit. Fig.|2]presents the best-fit composite spectra overplotted onto 
the observed ones, which were smoothed for clarity. The inferred 
velocities are listed in Table [3] The uncertainties correspond to 
the standard deviation of the residuals of the observed velocities 
compared to the best-fit model described in Sect. |4] 

3.2. Analysis Using Near-Infrared Atlas 

We also explored a different method of analysis by comparing 
our observations against near-IR spectra of early-type s tars of 
similar resolution and S/N. The work by iHanson et al.l (l2005h 
(hereafter, H05) provides a near-infrared atlas of Galactic hot 
stars that comprises of 37 well known, OB-type stars observed 
with VFT/ISAAC and Subaru/IRCS with a resolution of R ~ 
8 000 - 12 000 and a mean S/N ~ 100 - 300. The spectral types 
range from 03 to B3 and sample most luminosity classes. The 
advantage of using the atlas over the FASTWIND templates lies 
in the availability of extra diagnostic lines such as Civ at 2.08 
gm and the blend of the CNO complex at 2.115 gm with Hei 
at 2.112 jum. As mentioned in H05, the atlas should not be used 
for an accurate spectral classification but rather, to obtain a first 
estimate for the temperatures. 

Given that the H05 spectra are calibrated to air wavelengths 
(M. Hanson, private communication), we first applied an air-to- 
vacuum unit conversion. The resolution of the spectra was re¬ 
duced to match that of our data and the light ratio was fixed in 
accordance with the light curve analysis. We found that a com¬ 
posite spectrum of two identical 06.5 III spectra yielded an opti¬ 
mal fit to our observed spectra. However, there is no 06.5 main- 
sequence counterpart in the atlas, hence luminosity class at that 
temperature is not well-determined. The resulting spectral type 
is in agreement with the effective temperature of 38 kK adopted 
using FASTWIND models; the recent temperature-spectral type 


calibration for O-type stars bv ISimon-Dfaz et ^ d2014l) assigns 
an effective temperature of 38.7 kK to an 06.5 Milky Way dwarf 
of surface gravity 3.8 dex, with typical uncertainties of up to 1.5 
kK and 0.15 dex, respectively. 

The 06.5 III spectrum from H05 displays helium and hydro¬ 
gen features in both the H and K-band. While the centroid of the 
He I lines appear within 0.5 A from their theoretical wavelengths 
(in air, 1.70025 fim and 2.11201 /tm), the Hen line (2.18852 
jum) clearly displays a~ 1.5 A offset redward. Comparing to the 
adopted FASTWIND template, the blue wing of Hen appears 
weaker with respect to the red wing. It is unlikely that this is 
due to inefficient telluric correction, as the //-band He ii feature 
(1.691 84 fim) also suffers a similar, although smaller (~ 1 A) 
offset. In addition, a phase-dependent velocity discrepancy be¬ 
tween the CNO/He i blend and the He ii line is prominent in the 
ISAAC spectra, as can be seen in Fig. [3 The He ii line is mainly 
formed i n the transition region from the stellar photosphere to 
the wind (IFenorzer et al.ll2004ll . hence it is dominated by the un¬ 
certainty in the wind density structure. We therefore conclude 
that contamination by stellar winds is the most plausible expla¬ 
nation for this He ii offset (A. Herrero, private communication) 
as stellar winds are prominent in early-type stars. To avoid a def¬ 
inite bias of this effect to our study, we proceeded to measure 
velocities with the H05 spectra in two ways; first, fitting all fea¬ 
tures except for the Bry/He i blend and second, fitting only the 
2.112 fj.m He i line along with the CNO complex. While the first 
method relies on the contribution of four independent features, 
it is hindered by the uncertainty of the aforementioned He ii ef¬ 
fect and of a possible, inefficient correction in our observations 
for the strong, telluric CO 2 absorption band near 2.08 gm, in 
the region of C iv. Furthermore, both C iv and He ii are found at 
the edges of the observed regime where the S/N is expected to 
be lower and in addition, the blueshifted component of C iv is 
at most phases beyond the edge of our observed spectrum. Our 
second method employs fewer features, which, however, have a 
better S/N and do not suffer from telluric contamination. The de¬ 
rived velocities following both methods are listed in Table |4] as 
in Table [3] whereas the best-fit composite spectra are shown in 
the two panels of Fig. [3] For brevity, the method that excludes 
the Bry/Hei feature is labeled in the figures and tables as "- 
Bry/Hei", whereas the method that employs only the CNO/He i 
blend feature is labeled as "CNO/He i". 


4. Binary modeling 

We proceeded to a simultaneous fit of both the /c-band light 
curve and radial velocity curves with a detached, binary model. 
We used PHOEBE Subversion (release date, 2012-07-08) 
dPrsa & Zwitteill2005h. which implemen ts the Wilson-Devinney 
(WDl code (IWilson & DevinnevI [197 111 , to converge into the 
global optimization using Differential Corrections (DC) powered 
by a Fevenberg-Marquardt fitting scheme. We employed radial 
velocities measured using FASTWIND templates and the spectra 
from H05 following both methods discussed in Sec. 13.21 

The primary component is defined to be the star eclipsed at 
phase zero. The period determined from the AoV analysis of the 
photometric curve was assigned as an initial guess to the period 
in PHOEBE. We then fit the light curve with the following free 
parameters: time of primary eclipse HJDq, period P (days), in¬ 
clination i (deg), effective temperature of the secondary compo¬ 
nent Teff 2 (K) and surface potentials Oi, 2 - We fixed the effective 
temperature of the primary component at 38 kK. The light curve 
does not provide evidence for eccentricity, hence we fixed e - 0. 
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Both components were assumed to be tidally sychronized and 
surface albedos and gravity brightening exponents were fixed to 
unity, as for stars with ra diative envelopes. L imb darkening co¬ 
efficients were taken from Ivan Hammel (Il993li using the square- 
root law. The well-determined period and HJDq were fixed to 
their converged values and we fit the radial velocity curve for 
the semi-major axis a (R©), systemic velocity j (km s“^) and the 
ratio q of the mass of the secondary component over that of the 
primary compo nent. In addition, the Rossiter-McLaughlin effect 
(lRossiteijll924li was taken into account to correct for velocity 
shifts near the conjuction that occur when part of the approach¬ 
ing/receding surface of the occulted star is blocked. The root- 
mean-square (RMS) value of the fit of each velocity curve was 
assigned to the "sigma" parameter in PHOEBE, to weight the 
data in the overall cost function of the analysis. As a next step, 
both light and radial velocity curves were fit simultaneously, set¬ 
ting free all nine mentioned parameters. 

To prevent convergence to a local minimum of the solution 
space, we used the method of parameter-kicking. We first de¬ 
fined convergence as three consecutive iterations where adjusted 
parameters are less than or equal to their returned uncertainties. 
When the criterion is satisfied, all parameters are offset by a 
relative sigma ("kick"), defined as 


^^kick — 0.5 


100 


dPrsa & Zwitteil[2005l) . where Atot is the total number of points 
from the fitted photometric and velocity data-sets. Of the 1 000 
sets of parameters, we chose the set that yields the lowest sum 
values resulting from the fit of both the light and radial ve¬ 
locity curves. We then imported the particular set as input to a 
new run of 1 000 iterations, having the parameter-kicking option 
disabled. The mean and standard deviation of the values for each 
parameter were derived to provide the solution and the 1 -cr un¬ 
certainties, respectively. The above procedure was repeated three 
times, for each radial velocity set obtained with FASTWIND and 
the two methods for the 06.5 III spectrum from H05. The final 
values are presented in Table |5] It appears that the fit using He ii 
contributes to a separation of the components that is larger by 1 
Rq, compared to the method fitting only on the CNO/Hei fea¬ 
ture. In addition, complementing the He ii fit with that of C iv, 
yields components of almost equal mass. The systemic velocity 
modeled with the fit of the synthetic He ii from FASTWIND, is 
larger by ~ 20 km s ' than that modeled with the H05 spectra. 
The residuals of the three different fits are shown in Tables^and 
|4] We found that the velocity measurements based exclusively 
on the fit of the CNO/He i blend, yield the best-fit model as can 
be seen from the resulting RMS values of 11 and 8 km s * for 
the radial velocity curves of the primary and the secondary com¬ 
ponent, respectively. This best-fit model is displayed in Fig.|4] 
The physical parameters of D2-EB are presented in Table |6l 
All methods resulted in a surface gravity value for both com¬ 
ponents of log g ~ 3.9, which is typical for Galactic O-type 
dwarfs, whereas effective temperatures of 38 kK and 37 kK in¬ 
dicate spectral types of 06.5 and 07 respectively, with a typi- 
cal uncertainty of half a spectral type (iMartins et al.ll2005h . We 
adopted a conseryatiye uncertainty of IkK for the temperatures, 
equal to the step of the FASTWIND templates used. The inferred 
light ratio for each method was used to calculate the contribu¬ 
tion of the component spectra to the total composite spectrum, 
as mentioned in Section[3] In addition, we deriyed a negatiye fill¬ 
ing factor F -1.8 thus supporting the interpretation that both 

components are uneyolyed. The measured radii, period and in¬ 


clination yield synchronous rotational yelocities of ~ 140 and 
~ 130 km s“*, which are in good agreement with those inferred 
from the fitting process (Section [3TT l. 

Both methods that employ yelocities measured from the He ii 
line, yield a total mass which is ~ 10% greater than the method 
which discards He ii, owing to the larger separation of the com¬ 
ponents in the obseryed spectra. The latter best-fit method pro- 
yides measurements of the mass and radius of 24.42 + 0.15 M© 
and 9.69 + 0.07 R© for the primary and 22.06 + 0.68 M© and 
8.87 + 0.10 R© for the secondary component. The radii measure¬ 
ments correspond to radii of spheres with equal yolume to those 
defined by the WD surface potentials. Compared to the latest 
update d catalogue of well-studied, detached EBs by lSouthworthI 
(l2014h . the primary component of D2-EB is the third most mas- 
siye star to be studied with an precision better than ~ 2% and the 
most well-determined star with mass > 14.5 M©. Neyertheless, 
we caution that our 0.6% precision for the mass of the primary is 
deriyed as the uncertainty of the solution that optimizes the fit to 
our obseryations. Confirming this precision requires more yeloc- 
ity diagnostics and a larger number of obseryations taken at both 
quadratures. We hence conclude that the 0.6% precision should 
be taken with caution. For this reason, a different code than WD 
was also applied to eyaluate the current solution and reyise the 
uncertainties. 

We employed the genetic optimizer of EEC 
(lOrosz & Hauschiid^ l2000l) to fit the orbit based on eight 
input parameters. Specifically, the fit parameters in EEC were 
HJDq, P, i, the Roche lobe filling factor for each star, /, 
which is the ratio of the radius of the star toward the inner 
Lagrangian point Li to the distance to Li from the center 
of the star, / s jiCpointZ-^Li, the primary star’s yelocity 

semi-amplitude, Ki, and q. In fitting the orbit for D2-EB, EEC 
computed ~ 4 x 10'' orbits where the yalues of the eight input 
parameters yaried between fixed ranges judged to be applicable 
based on the WD fit of the CNO/He i radial yelocity set. The 
subsequent space was then projected as a function of each 
orbital and astro physical p a ramete r of interest in the same way 
as was done in IWilliamsI (l2009l) . From the global xtun^ 
estimate 1-cr uncertainties for deriyed and fitted parameters 
from the region where < xl„j„ + 1 ■ These yalues are listed in 
Table [6] The inferred masses and radii are 24.5 + 0.9 M© and 
9.2 + 0.1 R© for the primary and 21.7 + 0.8 M© and 8.7 + 0.1 R© 
for the secondary component. The best-fit model with EEC is 
displayed in Fig. |5] 

The masses deriyed from the best EEC fit are found to be 
consistent within uncertainties with those deriyed by PHOEBE 
and the WD code. The reyised uncertainties of the masses of the 
primary (3.7%) and the secondary (3.7%) are larger than those 
from our PHOEBE fit while the yalues of radii are lower by 2 - 
5% than those deriyed by the WD code and precise to 1%. This 
is not surprising, as the radii rely on the well-constrained light 
curye, while the masses rely more on the sparse radial yelocity 
data set. We adopt the yalues from the EEC analysis as their 
uncertainties are more conseryatiye. 


5. Distance 

A prerequisite to precisely determining independent distances 
to a double-lined eclipsing system is to fit a spectral energy 
distribution (SED) to accurate, multi-b and photometry taken 
at a known phase. iBaume et al.l (l2009ll conducted wide field 
UBVIc obseryations of the complex that hosts the Danks clusters 
and their surrounding field and obtained point spread function 
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(PSF) photometry of ~35 000 sources including our D2-EB0. 
The f/-band photometry was not available for D2-EB likely due 
to the high extinction. Using the ephemerides derived by our four 
methods, we c onverted the Julian d ates of the six per-filter expo¬ 
sures from the iBaume et al.l (l2009h data (G. Carraro, priv. com¬ 
munication) to units of phase. We sup plemented the op tical pho¬ 
tometry with 2MASS measurements dCutri ^ ^l2003h in the J, 
H, /fj-band obtained on JD 2 451 594.871^. All observations 
are consistent within their uncertainties, to out-of-eclipse phases 
of equal brightness (Eigs. SI IS. To convert t he ma gnitudes to 
fluxes, we used zeropoints from|Bessell et ^ (Il998h for the op¬ 
tical photometry and from ICohen et al.l (12003 ) for 2MASS. 

SEDs in accordance with the physical parameters of the best- 
fit spectra were generated with EASTWIND, to provide the flux 
density per surface unit through our studied bands. The compos¬ 
ite flux measured at Earth from a binary at a distance d, at a 
wavelength A, reddened to extinction A(/l), is given by 

fA = ^(R]Fi,,+rIF2,a) X 

where R 12 and Ui 2 are the radii and the surface fluxes of the 
components, respectively. The composite SED was reddened ac¬ 
cording to the new family of optic al and near-infrared ex t inctio n 
laws for O-type stars provided bv iMafz Aoellaniz et al.l (l2014ll . 
which c onstitute an improve ment of the widely used extinction 
laws bv ICardelli et al.l (Il989h . We ran a fitting algorithm over a 
wide range of distances with a step of 0.05 kpc, setting free the 
monochromatic parameters 7^5495 and £(4405-5495) for the type 
and amount of extinction respectively, and the best-fit values 
were considered to be those that minimized the weighted-;^^^. To 
estimate the uncertainty of our measurements, we used a Monte 
Carlo approach. In particular, we ran the fitting procedure 1 000 
times using sets of randomly selected parameters (photometry 
and radii) within their uncertainties assuming they are Gaussian- 
distributed. The corresponding values of distance are shown in 
Table | 6 ] for every set of radii and temperatures determined from 
the four different fit models. Our adopted radii measured with 
EEC yielded d = 3.52 + 0.08 kpc, £(4405 - 5495) = 3.66 + 0.06 
mag, £5495 = 3.26 + 0.04 and A 5495 = 11.9 + 0.1 mag, based 
on a photometric Te/y 2 ~ 36 kK. Assuming a spectroscopic 
Teff 2 - 37 kK, the distance changed slightly to r/ = 3.55 + 0.08 
kpc. Having three degrees of freedom, the reduced-;^f^ of our re¬ 
sulting fits was measured to be ~ 11 . 

We repeated the above procedure setting the amount of ex¬ 
tinction £(4405 - 5495) to the value of its band-integrated equiv¬ 
alent E(B - V) (and so increasing the degrees of freedom by 
one). Specifically, we calculated (B - y)o = -0.27 mag from the 
synthetic unreddened model at the isophotal wavelengths, which 
is in agree ment with the int r insic color of 06-9 giants/dwarfs 
provided by iMartins & PlezI (l2006h . thus yielding E(B -V)- 
3.28 + 0.16 mag. Our best-fit model then yielded d = 3.34 ±0.12 
kpc, though resulting from a less good fit (x^^^ ~ 23) than that 
achieved when setting £(4405-5495) free. Both reddened SEDs 
are shown in Eig. | 6 ] Adopting £(4405 - 5495) = 3.66 mag, 
the residuals indicate a better fit to the JHKs-hand photome¬ 
try and a reasonable fit to the U-band photometry. The £-band 
photometry clearly deviates from our adopted model causing 
the discrepancy of ~ 0.4 mag between the two values of the 
amount of extinction. This could be partly attributed to the un¬ 
certainty of the rather faint £-band photometric data (> 22 mag). 
However, it may also leave room for a further improvement of 

3 B = 22.41 ±0.16 mag, V = 19.40 ± 0.02 mag, £ = 14.73 ±0.03 mag 
11.42 ±0.03 mag,10.27 ± 0.03 mag, K, = 9.75 ± 0.03 mag 


the extinction laws. Indeed, having the B-band photometry ex¬ 
cluded from the fit, our procedure yielded d - 3.53 ± 0.08 and 
Ared ~ (with two degrees of freedom), which still deviates 
from;^^^^^ = 1. We caution that in the optical, b oth families of ex 


tinction laws by iMalz Anellaniz et alJ d2014t) and ICardelli et al.l 
(Il989h have been tested on a sample of low/intermediate red¬ 
dened stars (£(4405 - 5495) <1.5 mag). Revision of the optical 
laws for more optically-obscured targets, combined with a differ¬ 
ent value of the power law exponent for the near-infrared range, 
could result in a better fit. 

The measured distance to the D2-EB is in good agree¬ 
ment with reported values of di stance to the Dank s 2 cluster 
in the literature: 3.7 + 0.5 kpc (IChene et al. 2012h. 3 .8 ± 0.6 
kpc (iDavies et al.ll2012ll . 3.4 ± 0.2 knc dBica et akl^04ll . Never¬ 
theless, our inferred, well-constrained physical parameters pro¬ 
vide a precision of ~ 2%, which is a fa ctor of 3-8 improve - 
ment compared to the previous studies. iBaume et~^ (1200911 
showed that the cluster suffers from substantial reddening with 
E{B - V) - 2.4 mag (thus yielding Ay = 8.7 mag), which 
might be differential across their studied region. D2-EB resides 
outside their adopted boundary of Danks 2 and likely coincides 
with molecular gas that contaminates the north-east of the clus¬ 
ter (Eig. |7]), which appears to be the reason for our higher value 
of A 5495 = 11.9 mag. We further examined a Bpifzer/GLIMPSE 
map at 5.8 jum and found emission in the region of the gas. 


6. Discussion 

To measure the age of D2-EB and theref ore the Danks 2 clus- 
ter, we used the evolutionary models of lEkstrom et al.l (12012h 
for single stars, at Z = 0.014. The models assume initial stel¬ 
lar rotation at 40 % of the critical ve l ocity a ccounting for ef¬ 
fects discussed in iMevnet & Maede J (l2000l) . and correspond 
to main sequence velocities of 110 - 220 km s“', which are 
suitable for our case. All methods followed to derive masses 
agree with an age of ~ 3.2 Myr for both components, accord¬ 
ing to the Hertzsprung-Russell (H-R) diagram shown in Eig. [ 8 ] 
A by-eye interpolation of the evolutionary tracks implies that 
the members of D2-EB appear overluminous for their masses. 
In particular, using He ii for measuring velocities, the obtained 
masses found are ~ 15% lower than those predicted, for both 
components. When excluding Hen from the velocity analysis, 
this mass discrepancy increases to ~ 20%. A mass discrepancy 
is often reported when dealing with contact and semidetached 
massive binary systems jBurkholder et al.lll997t lBonanosll2009l : 

iKoumnia & Bonanosll2012h where one or both components are 

ovefllling their Roche lobe and is clearly attributed to mass trans¬ 
fer. Nevertheless, in the case of the well-detached D2-EB, we 
have no evidence for Roche lobe overflow and in addition, both 
components appear less massive than expected. Similar cases of 
ma ss discrepancy i n deta c hed systems have been also reported 
by iGonzalez et alJ ( 2005h. Williams et al.l (l2008h and most re- 
centlv bv iMassev et alJ (l2012l) . The latter study concluded that 
models with enhanced convective overshooting or higher initial 
rotation should be adopted in order to explain the ~ 11 % mass 
discrepancy observed in two det ached systems in the Lar ge Mag¬ 
ellanic Cloud. A recent study bv iMarkova & Puli (12015l) showed 
that a mass discrepancy for Galactic O stars with initial mass < 
35 Mq is evident, with the spectroscopically derived masses be¬ 
ing systematically lower than the evolutionary masses inferred 
from the theoretical models. Stars with spectroscopically mea¬ 
sured masses of ~ 25 Mq and ~ 22 Mq, equal to our adopted 
values for the components of D2-EB, were shown to display evo¬ 
lutionary masses of ~ 30 Mq and ~ 27 Mq respectively, which 
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are in a good agreement with our predictions from the theoretical 
mo dels. _ 

iDavies et alJ (1201 2h suggested an age for the cluster of 2 - 6 
Myr, based on the presence of the carbon-rich Wolf-Rayet D2- 
3 and the lumi n osity of the brightest OB supergiant in Danks 
2 . IChene et~^ (1201 2h fit theoretical isochrones to a colour- 
magnitude diagram of their near-infrared photometry and ar¬ 
gued that Danks 2 is likely older, with an age of 4 - 7 Myr. 
A temperature-independent method to evaluate the evolutionary 
status of D2-EB is the use of the mass-radius diagram, presented 
in Fig. |9] The two methods based on the fit of He ii are both con¬ 
sistent with an age of 4.5 Myr. Our adopted, best-fit model yields 
an age of ~ 5 Myr, following from both the WD code or ELC. 
All measured values of age are found to be within the range of 
ages determined for Danks 2 by the previous studies. 

D2-E B resides slightly outside the 1'5 + 0^5 angular size of 
Danks 2 (IChene et al.ll2012[ see Fig. |7j solid and dashed circles) 
beyond which, the cluster density profile falls below that of the 



cluster as a runaway binary. The radial velocity of D2-EB rela¬ 
tive to the cluster is 37 + 8 km s“'. This can be taken as a lower 
limit for the space velocity at which the system is escaping from 
the cluster. The angular distance 2'.1 from the center corresponds 
to a projected linear distance of 2.2 + 0.1 pc. Assuming a kine¬ 
matic age ~ 5 Myr equal to the measured age of the system, the 
tan gential velocity would be less than 1 km s '. 

iBlaauwl (1196 Ih hypothesized that runaway stars are created 
when a massive component within a binary system explodes as 
a Type II supernova. The secondary is then ejected with a ve¬ 
locity comparable to the orbital velocity at the time of t he su¬ 
perno va event. The explosion may not disrupt the system dHillsI 
Il983h . which will be observed as a OB-neutron star/black-hole 
system and eventually become a high-mass X-ray binary. Given 
that D2-EB is a double-lined O-type binary, it is unlikely that 
such a mechanism took place. Alternatively, D2-EB could have 
been part of a triple system including a very massive component 
(> 85 Mq) with a lifetime of less than 4.5 Myr. However, the 
inferred radial velocity curve lacks evidence of a low-mass third 
companion, unless it orbits with a high eccentricity and/or long 
period. Observations of D2-EB over a longer time span than the 
data presented here are req uired to invest i gate th e possibility of 
a third body in the system. iPoveda et al.l (Il967h suggested that 
runaway stars are dynamically ejected due to encounters of col¬ 
lapsing protostars in/near the core of young clusters. Interactions 
in clusters that contain initial "hard" binaries in crease the num - 
ber of escapees via binary-binary interactions (lMikkolalll983h . 
The less massive binary system reaches peculiar space velocities 
up to ~ 200 km s“*, while the more massive system trav els at 
a speed less than ~ 100 km s * ((Leonard & Duncanlll988h . The 
binary frequency for runaways with Voo > 30 km s ' is predicted 
to be 10% and it is striking that systems hosting components of 
M ~ 20 Mq as the most massive members are not predicte d 
to escape with more than 50 km s ' (iLeonard & Dunc^ll990h . 
This does not conflict with the observed lower limit for the space 
velocity of D2-EB. A young, dense cluster core increases the 
possibility of dynamical ejection, with the consequence that the 
kin ematic age of the runaw ays is similar to the age of the clus¬ 
ter dGualandris et al.l [20041) . In this case, the unreasonably low 
estimated value for the tangential velocity of D2-EB renders the 
possibility of the dynamical ejection through binary-binary in¬ 
teraction to be unlikely. 


7. Conclusions 

We present an analysis of new X-band spectra from VLT/ISAAC 
with which we determined accurate fundamental parameters of 
D2-EB, a massive, early-type eclipsing binary in the young clus¬ 
ter Danks 2, which is embedded in the G305 Galactic, star¬ 
forming region. The best-fit model to the binary was obtained by 
using the Hei line (2.112 pm) and the CNO complex for mea¬ 
suring radial velocities. The system was found to contain two 
co-evolutionary 06.5-7 main-sequence components with an age 
of ~ 5 Myr. We determined masses and radii of 24.5 ± 0.9 M© 
and 9.2 ± 0.1 R© for the primary and 21.7 + 0.8 M© and 8.7 ±0.1 
R© for the secondary component with a precision ~ 3.8% for the 
masses and ~ 1% for the radii. Models utilising a fit of Hen for 
measuring velocities, were found to yield a ~ 10% higher total 
mass and ~ 0.5 Myr younger age. 

Employing the precise measurements of the radii and out- 
of-eclipse optical and near-infrared photometry of D2-EB, we 
determined a distance to the system of c/ = 3.52 ± 0.08 kpc from 
a fit to the SED of the system. Up to now, this is the most well- 
constrained distance measurement to the Danks clusters and thus 
to the host complex, with a precision of ~ 2%. D2-EB is found 
to reside slightly outside the cluster and is optically obscured 
by molecular gas, which causes an extinction for the system of 
^5495 = 11.9 + 0.1 mag. We report a systemic velocity for the 
system that is inconsistent with that of Danks 2, making D2-EB 
a candidate runaway binary. However, neither of the two main 
mechanisms thought to yield runaways provide a robust expla¬ 
nation of the nature of the ejection of D2-EB from the Danks 2 
cluster. Future proper motion measurements with Gaifl which 
will be complete down to 20th mag in the visual, should con¬ 
strain the motion of both the cluster and the system, and provide 
well-determined values of the space velocity and the kinematic 
age of D2-EB. 
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Fig. 1. Reduced-^^ panel of the fit of composite FASTWIND spectra to the ISAAC spectrum of HJD 2456 345.856 71. Axes denote the effective 
temperature of the components and are split into six, equal intervals with 16 subdivisions of 10 km s“' each, to represent the projected rotational 
velocity range 50 - 200 km s“*. We find that the best fit corresponds to a pair of templates with Te/fi = 38 kK, Teff 2 = 37 kK and v sin i values of 
150 ± 20 and 130 ± 20 km s“', respectively. 



Wavelength (/rm) 

Fig. 2. Normalized ISAAC spectra (blue) sorted by the orbital phase (0) of D2-EB, including the Bry (2.166 yum)/Hei (2.162 pm) blend and the 
Hen (2.189 pm) line. The composite spectra (red) of two FASTWIND templates of 38kK and 37kK, weighted, summed and shifted to our derived 
velocities, are overplotted. The best-fit velocites are obtained from the fit of He ii. For clarity, a smoothing window is applied to the observed 
spectra and a constant is added to the relative fluxes. Absorption features in the wavelength range 2.170 - 2.185 pm are due to poor telluric 
correction. 
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2.115 


Fig. 3. Same as in Fig.[2l but the composite spectrum is comprised of two 06.5 III spectra taken from iHanson et al] (l200.5h . In the left panel, we 
show the entire wavelength range of the ISAAC spectra, which includes lines of Civ (2.08 /im), Hei (2.113, 2.162 yum), the CNO complex (2.115 
pm), Bry (2.166 pm) and Hen (2.189 pm). The best-fit velocities are measured by fitting all features except for the Bry/Hei blend. In the right 
panel, we show the result of the separate fit to the CNO/Hei blend, which was employed to derive velocities. 



Phase Phase 

Fig. 4. The left panel shows the /c-band light curve of D2-EB, phased to the best-fit parameters derived with PHOEBE using spectra from the H05 
atlas. Radial velocities were inferred from the fit of the CNO/He i blend. The right panel shows the radial velocity measurements overplotted by 
the modeled curve. The solid line denotes the primary and the dashed line the secondary component. Uncertainties correspond to the rms of the 
be st-fit, i.e. 11 km f or the primary and 8 km for the secondary. Green, solid lines represent the two observation dates of the BVf photometry 
bv iBaume et al.l ll2009h (six per filter), while the black dashed line indicates the date of the 2MASS measurements. Residuals of the fits are shown 
in the lower panels. 
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Fig. 5. Same as Fig. |4] but for the best-fit parameters inferred using ELC. Uncertainties for the radial velocity measurements correspond to 14 
km s“* for the primary and 10 km s“* for the secondary component. 
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logA (A) 

Fig. 6. Fit of the reddened, composite SEDs using FASTWIND models to the BVIcJHK, photometry, setting the amount of extinction £(4405 - 
5495) as a free parameter (red) and fixed to the value of the observed, band-integrated equivalent E{B — F) (blue). The former method provided 
the best-fit model yielding d - 3.52 ± 0.08 kpc, £(4405 - 5495) = 3.66 ± 0.06 mag, £5495 = 3.26 ± 0.04 and A 5495 = 11.9 ± 0.1 mag. 

Table 1. Summary of the parameters of the Danks 2 cluster and D2-EB. 


Parameter 

Danks 2 

D2-EB (this work) 

Total Mass (Mq) 

3 000 +800 ‘ 

46.2+1.2 

Age (Myr) 


4-72 

~5 

Distance (kpc) 

3.4 + 0.23 

3.52 + 0.08 

Visual Extinction (mag) 

8.74 

11.9 

References: * 

Davies et al. 

(I 2 OI 2 I), ^Chene et al.l(|2012b. 

teica et alJ (l2004l). iBaume et al.l (12009) 
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Fig. 7. Digitized Sky Survey (DSS-2) map of Danks 2 in the red filter . D2-EB (blu e arrow ) is located outside the l;5 radius of the cluster (red solid 
circle) with an uncertainty of ±015 (red dashed circles), as defined bv IChene et alJ (l2012h . 



log (T,^^/K) 


Fig. 8. H-R diagram showing the comparison of the parameters of D2-EB modeled with the WD code and our three, different sets of radial 
velocities obtained with: the fit of He it of FASTWIND templates (red circles), the fit of all Ai-band features of H05 spectra except for the Bry/He i 
blend (black triangles), the fit of the CNO/He i blend of the H05 spectra (green sq uares). The last set is also modeled with ELC and parameters 
are shown with blue asterisks. Evolutionary tracks and isochrones for single stars dEkstrom et aLllMTlI) at Z = 0.014 with rotation, are shown. 
Isochrones from left to right, correspond to 3.2, 4 and 5 Myr. All methods yield co-evolutionary components with an age of ~ 3.2 Myr that appear 
overluminous for their masses. 
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Fig- 9- Mass-Radius (M-R) diagram showing the comparison of the parameters of D2-EB with evolutionary tracks and isochrones for single stars 
dEkstrom et al.ll20T^ at Z = 0.014 with rotation. The symbols denoting each method are the same as in Eig.[^ Isochrones from the bottom up, 
correspond to 3.2, 4, 5 and 6.3 Myr. We find that both methods based on a Hen fit for measuring velocities yield co-evolutionary components of 
age 4.5 Myr. The best-fit model based on the CNO/He i blend fit yields an age of 5 Myr, as modeled both with the WD code and EEC. 


Table 2. Log of spectroscopic observations of D2-EB with VLT/ISAAC. 


HJD 

Total integration time (s) 

Airmass 

2456 340.78201 

2 x100 

1.31 

2456 345.85671 

2 x100 

1.28 

2456 360.82868 

2 x100 

1.30 

2456 360.84070 

2 x100 

1.31 

2456 360.85982 

2 x100 

1.35 

2456 378.695 33 

2 x100 

1.29 

2456 378.80874* 

2 x100 

1.34 

2456 383.71590 

2 x100 

1.27 


(*) Discarded from our analysis. 


Table 3. Radial velocity measurements of D2-EB using EASTWIND models. 


HJD 

Orbital 

phase 

RVi 

(kms-‘) 

FASTWIND 
(O-C)i 
(km s“*) 

RV2 

(km s“*) 

(0-C)2 
(km s“*) 

2456 340.782 01 

0.19 

-181±19 

24 

279 ± 30 

32 

2456345.85671 

0.69 

229 ± 19 

-1 

-171 ±30 

55 

2456360.828 68 

0.13 

-173±19 

-12 

207 ± 30 

8 

2456 360.84070 

0.14 

-183±19 

-19 

182 ±30 

-21 

2456 360.859 82 

0.14 

-193±19 

-23 

217 ±30 

8 

2456378.695 33 

0.43 

-117± 19 

-27 

88 ±30 

-34 

2456 383.71590 

0.92 

142 ± 19 

13 

-98 ± 30 

18 
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Table 4. Radial velocity measurements of D2-EB using spectra from H05. 


HJD 

Orbital 

phase 

RVi 

(km s“*) 

H05 (-Bry/Hei) 
(O-C)i 
(km s *) 

RV2 

(km s“*) 

(0-C)2 
(km s *) 

Orbital 

phase 

RVi 

(km s“') 

H05 (CNO/Hei) 
(O-C)i 
(km s ') 

RV2 

(km s“*) 

(0-C)2 
(km s *) 

2456 340.782 01 

0.20 

-241±14 

1 

209 ± 11 

-12 

0.21 

-241±11 

-17 

229 ±8 

-4 

2456345.85671 

0.70 

189 ± 14 

-29 - 

241 ± 11 

8 

0.71 

194 ± 11 

-17 

-236 ± 8 

12 

2456 360.82868 

0.14 

-213±14 

-11 

197 ± 11 

17 

0.15 

-203 ±11 

-13 

207 ±8 

11 

2456 360.84070 

0.15 

-223 ±14 

-18 

197 ± 11 

13 

0.15 

-188± 11 

5 

207 ±8 

8 

2456360.859 82 

0.15 

-213±14 

-3 

182 ± 11 

-7 

0.16 

-203 ±11 

-6 

197 ±8 

-7 

2456 378.695 33 

0.44 

-117± 14 

-14 

83± 11 

3 

0.45 

-87± 11 

-8 

73 ±8 

-2 

2456 383.715 90 

0.93 

107 ± 14 

14 - 

113 ± 11 

8 

0.94 

92± 11 

10 

-103 ± 8 

1 


Table 5. Results from the analysis of the light and radial velocity curves with PHOEBE. 


Parameter 

FASTWIND 

H05 (-Bry/He i) 

H05 (CNO/Hei) 

Period, P (days) 

3.373 652 ± 0.000 002 

3.373 474 ± 0.000 003 

3.373 335 ± 0.000003 

Time of primary eclipse, HJDq 

2 455 685.653 33 ± 0.0001 

2455 685.652 73 ± 0.0001 

2455 685.652 45 ±0.0001 

Inclination, i (deg) 

68.4 ± 0.2 

68.0 ±0.2 

68.2 ±0.2 

Surface potential, Oj 

4.50 ± 0.02 

4.63 ± 0.05 

4.47 ± 0.03 

Surface potential, 02 

4.62 ± 0.04 

4.68 ± 0.05 

4.56 ± 0.08 

Light ratio in I, L 2 /L 1 

0.79 ±0.01 

0.89 ± 0.05 

0.79 ± 0.01 

Mass ratio, q 

0.92 ± 0.01 

0.98 ± 0.02 

0.90 ± 0.02 

Systemic velocity, y (km s“^) 

11 ± 1 

-13± 1 

-7± 1 

Semi-major axis, a (R©) 

34.95 ± 0.04 

35.08 ±0.17 

34.02 ± 0.20 

Semi-amplitude, K\ (km s“^) 

233 ± 1 

241 ±4 

225 ±5 

Semi-amplitude, K 2 (km s“^) 

254 ± 1 

247 ±2 

249 ±3 

Radius, Ti^poie 

0.276 ± 0.002 

0.271 ±0.002 

0.278 ± 0.002 

. ^l,side 

0.282 ± 0.002 

0.277 ± 0.002 

0.284 ± 0.002 

. point 

0.295 ± 0.003 

0.290 ± 0.003 

0.297 ± 0.003 

. ^l,back 

0.290 ± 0.003 

0.285 ± 0.002 

0.292 ± 0.002 

Radius, r 2 ,poie 

0.254 ± 0.002 

0.264 ± 0.002 

0.255 ± 0.002 

. ^2,side 

0.259 ± 0.003 

0.269 ± 0.002 

0.260 ± 0.002 

. ^2,point 

0.269 ± 0.003 

0.280 ± 0.003 

0.271 ±0.003 

. ^2,back 

0.266 ± 0.003 

0.276 ± 0.003 

0.267 ± 0.003 


Table 6. Physical parameters of D2-EB. 


Parameter 

FASTWIND 

H05 (-Bry/Hei) 

H05 (CNO/He i) H05 (CNO/He i) (ELC) 

Ml (Mo) 

26.28 + 0.20 

25.80 + 0.09 

24.42 + 0.15 

24.5 + 0.9 

M2 (Mo) 

24.12 + 0.16 

25.16 + 0.66 

22.06 + 0.68 

21.7 + 0.8 

Ri (Ro) 

9.90 + 0.05 

9.75 + 0.07 

9.69 + 0.07 

9.2+ 0.1 

Ri (Ro) 

9.08 + 0.06 

9.46 + 0.19 

8.87 + 0.10 

8.7+ 0.1 

loggi 

3.87 + 0.01 

3.87 + 0.01 

3.85 + 0.01 

3.90 + 0.01 

log g 2 

3.90 + 0.01 

3.89 + 0.01 

3.89 + 0.01 

3.90 + 0.01 

Teffl (K) 


38000 (fixed) 


T.//2 (K) 

36833 + 1000 

36964 + 1000 

36788 + 1000 

36400 + 1000 

logLi/Lo 

5.26 + 0.05 

5.25 + 0.05 

5.24 + 0.05 

5.20 + 0.05 

log L 2 ILQ 

5.13 + 0.05 

5.18 + 0.05 

5.11+0.05 

5.08 + 0.05 

Distance (kpc) 

3.77 + 0.08 

3.81+0.09 

3.68 + 0.08 

3.52 + 0.08 
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